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ABSTRACT 

Laser  induced  evaporation  of  material  from  the  surface 
of  an  aluminum  target  in  a  vacuum  was  studied.   Based  on  a 
literature  examination,  material  removal  using  a  normal 
pulse  laser  was  judged  to  be  more  efficient  than  for  a 
Q-switched  laser.   The  experiment  was  conducted  using  a 

neodymium  glass  laser  modified  for  normal  pulse  operation. 
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The  energy  density  was  varied  from  8.5x10   J/cm  where  no 
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breakdown  occurred  to  5x10   J/cm  where  the  threshold  for 

breakdown  was  exceeded.   The  normal  pulse  duration  was 

600  jjs.   Analysis  of  the  ejected  material  was  achieved  by 

using  a  Hughes  Ionization  Gauge  placed  in  the  path  of  the 

ejected  material.   Oscilloscope  traces  of  the  ionization 

gauge  output  show  that  the  gauge  "sees"  what  is  flying  past 

it.   There  is  good  correlation  between  laser  radiation, 

plasma  radiation  and  ionization  gauge  fluctuations.   The 

ionization  gauge  gave  distinguishable  signals  for  ions, 

electrons,  and  neutral  particles  ejected  from  the  target 

surface.   Signal  sequence  was  dependent  on  the  particle 

velocity.   By  measuring  the  elapsed  time  after  ejection 

from  the  surface  and  the  target  to  collector  distance,  the 

first  arriving  neutral  particle  velocity  was  determined  to 
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be  5.2x10   cm/s. 
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I.   INTRODUCTION 

A  wire  explosion  is  a  complex  phenomenon  resulting  from 
very  rapid  electrical  heating  of  a  piece  of  metal  which  has 
a  small  cross  sectional  area.   The  wire  is  heated  and  under- 
goes phase  change.   There  is  a  bright  flash  of  light  and  a 
loud  report  (in  air).   Thus,  the  term  "exploding  wire"  has 
been  coined  to  describe  this  phenomenon.   More  appropriately, 
the  term  "exploding  conductor"  will  be  used  to  denote  the 
explosion  of  a  single  wire,  a  multiple  wire  array,  a  metal 
foil,  or  a  gas  puff  by  a  high  current  pulse.   This  technique 
is  a  relatively  simple  way  of  generating  a  plasma,  which 
can  be  used  to  study  the  dense  plasma  itself  or  for  further 
experimental  work  such  as  spectroscopy. 

The  idea  of  exploding  conductors  is  not  a  new  one  [3] . 
As  early  as  1773,  Edward  Nairne  conducted  experiments 
exploding  a  .15  mils  diameter  iron  wire  in  proving  that 
current  in  all  parts  of  a  series  circuit  is  the  same  [45] . 
Eighty  years  later,  Michael  Faraday  reported  on  wire  explo- 
sions used  to  produce  a  metal  film  or  mirror.   The  use  of 
exploding  wires  was  of  little  importance  until  the  work  of 
John  A.  Anderson  in  the  1920 's.   At  Mt.  Wilson  Observatory, 
he  showed  that  the  temperatures  involved  approached  those 
of  the  sun,  in  excess  of  3000°C  [48] .   In  the  years  that 
followed,  research  on  exploding  conductors  was  principally 
conducted  by  Nabaoka,  Kleen,  Wrana,  Eiselt,  Conn,  Kvartskhava 


and  Lebedev  [3].   Since  1950,  exploding  conductors  have 
become  a  matter  of  great  scientific  interest.   Perhaps  the 
largest  number  of  researchers  to  use  exploding  conductors 
have  capitalized  on  the  extremely  short,  intense  light  out- 
put.  These  intense  pulses  of  light  have  been  used  in  spec- 
troscopy and  high  speed  photography,  as  well  as  in  radiation 
chemistry  and  laser  technology  as  an  excitation  source  for 
coherent  radiation  [50] . 

Other  uses  of  exploding  conductors  are  the  production 
and  study  of  aerosols,  production  of  thin  metal  films  [27], 
silvering  of  mirrors,  and  chemical  synthesis  [1] .   In  the 
past  10-12  years,  the  use  of  exploding  conductors  as  bridge- 
type  electrical  detonators  has  found  widespread  application 
both  as  a  method  for  multi-point  detonation  of  explosives 
[46]  and  as  a  method  to  generate  shock  waves  in  different 
materials  [9] .   Another  area  of  intense  research,  and  one 
in  which  this  paper  is  related  to,  is  the  use  of  exploding 
conductors  in  simulating  nuclear  weapons  effects.   In  recent 
years,  progress  has  been  made  in  the  development  of  exploding 
conductors  as  x-ray  (photon)  sources  for  use  in  simulating 
nuclear  weapons  effects  on  various  satellite  "black  boxes". 

The  emission  of  electromagnetic  radiation  is  a  result 
of  suddenly  electrically  heating  thin  wires.   Four  confer- 
ences held  between  1958  and  1967  dealt  exclusively  with 
this  subject  [17,18,19,20].   Toward  this  end,  several  varia- 
tions of  the  exploding  conductor  approach  to  high  energy 
photon  generation  have  been  explored  in  the  past  few  years. 


Single  and  multiple  wire  arrays  [8,47],  cylindrical  foils 
[60],  and  gas  puffs  have  all  been  tried  to  date  with  limited 
success.   The  energy  of  the  emitted  x-rays  has  been  too  low 
to  be  used  in  nuclear  explosion  simulations.   The  problem 
centers  around  insufficient  energy  being  coupled  to  the 
plasma  to  allow  for  radiation  in  the  x-ray  region. 

Whatever  the  application  or  shapes  of  the  conductor, 
the  electrical  circuit  and  the  components  are  similar  and 
perform  similar  functions,  the  only  difference  being  in 
size  and  characteristics.   Each  circuit  contains  an  energy 
storage  device,  a  switch  or  triggering  device,  and  an  explo- 
sive conductor.   The  energy  storage  device  may  be  a  bank 
of  large  capacitors  or  special  multiterawatt  generators, 
such  as  Pithon  at  Physics  International,  Gamble  II  at  the 
Naval  Research  Laboratory,  or  Blackjack  IV  at  Maxwell 
Laboratories,  Inc.,  [63].   These  devices  are  capable  of 
producing  100  KJ  pulses  of  50-75  nanosecond  duration.   This 
high  rate  of  energy  release  requires  special  switching 
devices  and  transmission  lines.   Basically,  a  Marx-capacitor 
bank  (200-300  KJ)  is  used  to  charge  a  single,  large,  water 
capacitor.   This  capacitor  discharges  into  a  water  dielec- 
tric coaxial  transmission  line  that  tapers  to  a  l-2fi  output 
resistance  at  the  point  where  the  pulse  is  extracted  into 
a  vacuum  diode.   Thin  wires,  ribbons  or  foils  are  stretched 
between  the  electrode  gap  in  the  diode  [41] .   In  the  case 
of  a  wire  (13,29,52,53,68]  or  a  foil  [66],  the  initial  part 


of  the  electrical  discharge  causes  the  conductor  to  rapidly 
vaporize  uniformly  along  its  entire  length.   The  plasma  so 
formed  is  then  compressed  in  a  Z-pinch.   As  a  result,  the 
plasma  consists  of  "pinched"  and  "flared"  structures  char- 
acteristic of  the  m  =  0  sausage  instability  to  which  Z- 
pinches  are  susceptible.   It  is  from  these  pinched  and  flared 
areas  that  radiation  is  emitted  [22].   Gas  puff  experiments 
conducted  at  Physics  International  are  similar  to  the 
exploding  conductor  ones  just  described  except  that  a  gas 
puff  is  substituted  for  the  wires.   In  addition,  the  gas 
puff  is  ionized  with  microwaves  prior  to  the  arrival  of  the 
electrical  pulse  to  insure  conductivity  of  the  gas. 

Most  of  the  methods  giving  rise  to  electromagnetic 
radiation  emitted  by  plasmas  are  atomic  in  nature,  that  is, 
they  are  due  to  transitions  involving  bound  or  free  states 
of  atoms  or  ions.   There  are  three  basic  methods  of  energy 
radiation:   bremsstrahlung,  recombination  radiation,  and 
discrete  radiation.   Bremsstrahlung  or  free-free  radiation 
occurs  when  a  free  electron  collides  with  an  ion  or  neutral 
particle  and  makes  a  transition  to  another  free  state  of 
lower  energy  with  the  emission  of  a  photon.   The  energy  of 
the  photon  may  amount  to  any  fraction  of  the  initial  kinetic 
energy  of  the  electron,  KT  ;  therefore,  the  bremsstrahlung 

has  a  continuum  of  frequencies  determined  by  KT  .   At 

4 
T  =  10   °K,  the  bremsstrahlung  lies  in  the  visible  and 

g 

infrared  regions,  while  at  T  =10   °K  most  of  the  radiation 
lies  in  the  x-ray  region.   Recombination  occurs  when  a 
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free  electron  is  captured  by  an  n-times  ionized  atom  and 
makes  a  transition  to  a  bound  state  of  the  (n-1)  times 
ionized  atom.   The  surplus  energy  is  emitted  as  recombina- 
tion or  free-bound  radiation.   The  energy  is  equal  to  the 
sum  of  the  kinetic  energy  of  the  free  electron,  KT  ,  and 
its  binding  energy.   Since  free  electrons  have  a  continuous 
energy  spectrum,  the  photons  emitted  form  a  continuous 
energy  spectrum.   Discrete  radiation  or  bound-bound  radiation 
occurs  when  a  bound  electron  of  an  atom  or  ion  is  excited 
to  a  higher  energy  state  by  particle  collisions.   When  the 
electron  returns  to  a  lower  state,  it  yields  energy  at  various 
discrete  frequencies.   Whenever  a  plasma  contains  ions  that 
have  not  been  completely  stripped  of  orbital  electrons,  this 
radiation  will  appear.   Generally,  it  takes  greater  plasma 
energy  to  completely  strip  atoms  or  molecules  with  higher 
atomic  numbers.   As  the  electron  temperature  is  increased 
more  tightly  bound  electrons  are  removed.   This  will  be 
associated  with  an  increase  in  the  excitation  energy  of  the 
ions.   As  a  result,  the  line  spectrum  will  shift  from  the 
visible  to  the  x-ray  region.   Although  the  radiation  emitted 
from  exploding  conductors  is  a  combination  of  the  above 
processes,  it  is  dominated  at  high  temperatures  by  bremsstrah- 
lung  and  discrete  radiation  in  the  x-ray  region  [30,34]. 

The  spectral  yields  and  hardness  of  the  emitted  radiation 
from  exploding  conductors  depends  on  the  ionic  state  and  the 
electron  temperature  which  in  turn  are  influenced  by  a 
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combination  of  competing  physical  processes.   Principal 
among  these  are  the  available  generator  power,  the  magneto- 
hydrodynamic  fluid  motion  of  the  wire  plasma,  the  coupling 
of  the  wire  plasma  to  the  driving  electrical  circuit  and  the 
atomic  physics  and  energetics  of  the  plasma.   The  z-number, 
geometry  and  mass  of  the  wire  material  also  play  an  impor- 
tant role  in  the  above  process.   A  good  deal  of  the  under- 
standing of  how  these  competing  processes  affect  the  dynamic 
and  radiative  behavior  of  the  wire  plasma  is  obtained  from 
experimental  data.   In  addition,  computer  codes  such  as 
WHYRAC  [51]  have  been  developed  to  calculate  the  detailed 
time  history  of  the  wire  implosion.   The  codes  provide  in- 
sight into  the  coupling  of  plasma  to  the  generators  and  into 
the  energetics  of  the  implosion.   Although  many  of  the 
processes  are  fairly  well  understood,  their  relationship 
and  affect  on  one  another  is  not. 

It  is  theorized  that  one  way  to  obtain  x  rays  with  the 
desired  characteristics  may  be  to  use  materials  with  differ- 
ent, higher  z  values  or  even  combinations  of  materials. 
Many  materials  can  not  be  fabricated  into  thin  wires  or 
do  not  have  a  high  enough  vapor  pressure  to  allow  the  forma- 
tion of  a  vapor.   In  response,  investigators  have  therefore 
proposed  that  the  vapor  puffs  be  generated  by  using  a  laser. 
One  of  the  vacuum  diodes  would  be  fabricated  out  of  the 
target  material.   The  target  electrode,  irradiated  by  the 
laser,  would  experience  a  steady  temperature  rise.   If  the 
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laser  beam  is  intense  enough,  the  surface  of  the  target 
would  melt  and  begin  to  vaporize.   The  vapor  expands, 
exposing  fresh  target  layers  to  heating  and  vaporization. 
Depending  on  the  laser  intensity,  the  vapor  may  remain 
un- ionized  or  a  plasma  may  be  generated.   The  high  voltage, 
high  current  pulse  from  the  Marx-generator  would  then  be 
timed  to  pass  through  the  cloud  of  plasma,  neutral  particles 
and  molten  metal.   This  procedure  advances  the  possibility 
of  using  virtually  any  material  or  combination  of  materials, 
in  addition  to  being  able  to  vary  the  mass  and  temperature 
of  the  vapor  cloud  by  varing  the  power  and  intensity  of  the 
laser  pulse. 

As  indicated  earlier,  one  important  parameter  on  which 
the  dynamics  of  the  z-pinch  and  x-ray  emission  depends  is 
the  mass  of  the  pinched  column.   This  thesis  investigates 
the  mass  of  material  which  is  blown  off  a  surface  as  a  func- 
tion of  the  laser  parameters  and  a  method  by  which  it  can 
be  measured.   Aluminum  was  selected  as  the  initial  target 
material  because  of  the  extensive  work  done  with  aluminum 
exploding  wires,  thus  facilitating  a  comparison  of  the  two 
methods  in  the  future . 
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II.   THEORY 

The  effect  caused  by  the  absorption  of  high  power  laser 
radiation  on  the  surface  of  an  opaque  solid  slab  includes 
surface  heating,  melting  and  vaporization.   If  the  laser 
intensity  is  low  enough  the  vapor  is  un-ionized  and  remains 
transparent  to  the  incoming  laser  radiation.   At  higher 
intensities,  the  vapor  becomes  ionized  generating  a  plasma. 
Throughout  this  thesis,  I  will  use  the  term  "vapor"  to  de- 
note an  ejected  cloud  of  material  which  contains  both  neu- 
tral particles  and  plasma.   In  conjunction  with  the  radiation 
impinging  on  the  opaque  surface,  a  cycle  of  events  occurs. 
This  is  depicted  in  Figure  1.   Prior  to  considering  these 
events  individually,  a  few  comments  on  lasers  in  general 
are  appropriate. 

First,  lasers  can  be  divided  into  two  categories,  pulsed 
and  continuous.   Only  pulsed  lasers  were  considered  in  this 
study.   Pulsed  lasers,  in  turn,  may  be  subdivided  into 
normal  pulse  and  Q-switched  depending  on  pulse  operation. 
A  normal  pulse  laser  is  one  which  is  pumped  by  a  flashlamp 
and  the  radiation  allowed  to  emerge  when  the  threshold  con- 
ditions for  lasing  are  reached.   Normal  pulse  lasers  emit 

in  pulses  lasting  tens  of  milliseconds  with  peak  powers  of 

4       5 
10   to  10  watts.   The  term  Q-switched  laser  denotes  lasers 

employing  an  element  of  variable  loss  within  the  cavity  and 

g 

emitting  peak  powers  greater  than  10   watts  and  lasting 
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Figure  1.  Events  diagram  for  laser-target  interaction.  Dotted  boxes 
are  considered  negligible  or  do  not  contribute  to  the  study 
at  hand  and  will  not  be  considered.  For  simplicity, 
feedback  mechanisms  between  processes  are  not  included. 
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several  tens  of  nanoseconds.   With  Q-switched  lasers,  flux 

9         2 
densities  in  excess  of  10   watts/cm   are  easily  attained. 

In  this  regime,  the  vaporization  temperature  of  any  metal 
will  be  reached  in  less  than  one  nanosecond.   At  this  point, 
the  input  energy  begins  to  supply  the  latent  heat  of  vapori- 
zation to  a  thin  layer  of  material  at  the  surface  and  break- 
down (avalanche  ionization)  will  occur.   During  breakdown, 
the  few  naturally  occurring  thermal  electrons  from  the  target 
surface  are  heated  by  the  laser  beam.   They  will  gain  energy 
until  they  are  capable  of  ionizing  neutral  particles  by 
collision.   Each  collision  will  produce  ions  and  more  elec- 
trons which  will  also  be  heated  by  the  laser  beam  to  the  point 
where  they  too  can  ionize  neutral  particles.   The  process 
thus  continues  exponentially  producing  electrons  and  ions 
and  absorbing  more  and  more  of  the  laser  radiation.   The 
"breakdown  threshold"  is  the  laser  intensity  required  to 
generate  a  high  density  plasma  at  the  target  surface.   Experi- 
ments have  shown  that  breakdown  over  a  laser- irradiated  sur- 
face exhibits  a  sharp  threshold  in  laser  intensity.   Several 
models  have  been  developed  that  predict  this  breakdown  thres- 
hold [40,54].   Below  threshold,  insufficient  energy  is  avail- 
able to  propogate  the  avalanche  ionization  and  the  vapor 
remains  relatively  transparent.   When  the  breakdown  threshold 
is  reached,  a  high  absorption  plasma  cloud  is  formed  on  the 
target  surface.   For  example, consider  Q-switched  neodymium- 
glass  laser  radiation  (1.06  ym)  incident  upon  an  aluminum 
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9       10 
target.   A  typical  pulse  could  have  a  flux  of  10   to  10 

2 
watts/cm   and  a  duration  of  10-100  nanoseconds.   Figure  2, 

taken  from  Ref.  4  models  the  relationship  between  the  thres- 
hold intensity  required  for  breakdown  and  the  laser  pulse 
duration.   The  prediction  is  that  breakdown  will  occur  in 
this  case.   Because  of  breakdown,  the  energy  in  the  beam  is 
devoted  to  heating  a  small  amount  of  vaporized  material  to 
a  high  temperature,  while  the  heat  transfer  from  the  hot 
vapor  to  the  bulk  of  the  target  material  is  limited  due  to 
the  short  pulse  duration.   A  normal  pulse  laser  has  a  series 
of  spikes,  many  of  which  are  capable  of  producing  breakdown. 
The  difference  is  that  the  individual  spikes  do  not  possess 
enough  intensity  to  perpetuate  the  ionization  process.   The 
vapor  is  heated,  breakdown  occurs  and  absoprtion  of  energy 
increases,  but  the  laser  spike  intensity  falls  off  rapidly 
and  breakdown  ceases.   This  process  is  repeated  many  times 
resulting  in  more  energy  being  coupled  to  the  target  and  thus 
more  material  being  ejected.   Consequently,  a  given  amount  of 
energy  delivered  at  very  high  power  is  less  effective  in 
causing  vaporization  than  the  same  amount  of  energy  delivered 
in  a  longer,  lower  power  pulse.   It  is  for  this  reason  that 
only  normal  pulse  lasers  were  considered  in  this  thesis. 

A.   REFLECTIVITY 

When  one  considers  the  coupling  of  laser  radiation  to 
the  surface  of  a  target,  one  first  needs  to  know  how  much 
energy  is  absorbed  into  the  material.   The  Drude-Lorentz  free 
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Figure  2.  Pulse  duration  time  vs  breakdown  threshold  intensity  for  an 
aluminum  target  irradiated  by  a  1.06  ym  laser  pulse. 
Alternatively,  this  graph  gives  the  breakdown  time  vs  peak 
value  of  incident  laser  beam  intensity.  Taken  from  Ref .  4. 
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electron  theory  provides  an  acceptable  model  for  metals 
interacting  with  infrared  wavelength  photons  [66].   According 
to  the  model,  electromagnetic  radiation  interacts  only  with 
free  electrons  in  a  metal.   The  absorbed  photon  energy  raises 
electrons  to  higher  energy  states  in  the  conduction  band. 
The  excited  electrons,  in  turn,  collide  with  lattice  phonons 
and  other  electrons  as  they  give  up  their  energy.   The  electron- 
phonon  collision  frequency  is  proportional  to  the  phonon  popu- 
lation in  the  metal.   The  phonon  population,  on  the  other  hand, 
determines  the  temperature  of  the  metal.   Thus,  the  rise  in 
temperature  strongly  affects  the  electron-phonon  collision 
frequency  which  in  turn  affects  the  reflectivity  of  the  metal. 
Based  on  the  Drude-Lorentz  theory  several  authors  have  developed 
theoretical  expressions  that  predict  the  absorbtivity  of 
metals  [49,66] . 

The  reflectivity  of  a  "real  metal"  surface  is  largely  an 
empirical  matter.   Bonch-Bruevich  et  al.  [10]  investigated 
the  reflectivity  at  1.06  urn  of  aluminum,  copper,  dural,  steel, 
and  silver  as  the  metals  were  irradiated  with  a  laser  beam. 
The  investigators  surrounded  the  sample  with  a  sphere  to 
monitor  the  reflected  radiation,  Figure  3.   Pronounced  changes 
both  during  each  individual  spike  and  over  the  pulse  as  a 
whole  were  reported.   Figure  4  shows  the  generalized  reflec- 
tivity behavior  of  a  metal  during  an  individual  spike.   The 
region  of  rapid  change  of  reflectivity  a-b  is  due  to  initial  heating 
of  the  metal.   The  segment  b-c  has  a  constant  reflectivity  which  the 
authors  attributed  to  the  constant  temperature  during  the  time 
when  the  melting  wave  is  propagating  into  the  metal.   Zavecz  [71], 
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Figure  3.   Schematic  of  the  apparatus  used  by  Bonch-Bruevich, 
et  al.  [10]  for  measuring  the  change  in  reflecting 
power  of  a  metal  under  the  action  of  laser 
radiation . 
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Figure  4.   Change  of  reflecting  power  of  silver  at  1.06  ym 

during  a  single  spike.   The  form  of  the  laser  spike 
with  energy  of  7.5  KJ/cm2  is  also  depicted  [10]. 
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disagreeing,  feels  that  this  low  reflectivity  plateau  is 
a  function  of  the  instantaneous  pulse  intensity  and  the 
specimen's  equilibrium  temperature.   Finally,  the  surface 
temperature  begins  to  rise  again,  thus  causing  a  further  drop 
in  reflectivity,  segment  c-d.   The  increase  in  reflectivity 
to  the  right  of  point  d  is  due  to  the  decay  of  the  spike 
flux  density  and  the  subsequent  drop  in  temperature.   Similar 
results  were  published  by  Chun  and  Rose  [21]  as  reflected  in 
Figure  5.   In  addition,  Chun  and  Rose  examined  the  dependence 
of  absorptance  on  the  depth  (internal  volume)  of  the  laser 
crater.   These  results,  depicted  in  Figure  6,  indicate  that 
absorptance  changes  as  the  depth  of  the  crater  changes.   This 
means  that  even  if  the  initial  reflectivity  is  high,  much  of 
the  energy  in  a  focused  normal  pulse  laser  will  be  absorbed 
by  a  metal  surface.   Although  only  experimental  results  for 
incident  radiation  at  1.06  um  has  been  described,  experimental 
work  at  other  wavelengths  has  been  done  [57] .   In  all  cases 
the  experimental  observations  show  that  under  high-intensity 
radiation,  metal  absorptivity  increases  non  linearly. 

B.   HEATING  OF  THE  MATERIAL 

Possibly  one  of  the  most  important  effects  of  intense 
laser  irradiation  is  the  conversion  of  optical  energy  in  the 
beam  into  thermal  energy  in  the  material.   The  radiation  mean 

free  path  for  visible  and  infrared  wavelengths  in  a  solid 

-4 
material  such  as  a  metal,  is  typically  of  the  order  of  10    cm 

or  less,  so  that  the  deposition  of  laser  energy  can  be  considered 
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Figure  5.   Time  dependent  reflectance  measurement  at 

1.06  ym  from  Ref.  21.   A  single  laser  pulse 
with  a  power  density  of  107  W/cm2  was  used. 
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Figure  6.   Dependence  of  absorbtance  at  1.06  um  on  crater 
volume  from  Ref.  21.   m/p  is  the  ratio  of 
mass  removed  to  material  density. 
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to  be  a  surface  phenomenon  insofar  as  the  transport  of  energy 
in  a  solid  is  concerned.   The  Drude-Lorentz  free  electron 
theory  as  previously  discussed  is  the  basis  of  thermal 
response  of  the  material.   The  time  required  (relaxation  time) 
for  electrons,  excited  by  irradiation,  to  transfer  their 
energy  to  the  lattice  by  means  of  electron-phonon  collisions 
amounts  to  approximately  10     sec  [36,72,74].   When  compared 
to  a  laser  pulse,  the  relaxation  time  is  short,  allowing  one 
to  assume  that  the  heat  transfer  to  the  solid  is  instantaneous 
because  local  equilibrium  is  rapidly  established.   Therefore, 
one  is  justified  in  assuming  that  temperature  is  a  valid 
concept  and  so  the  normal  equations  of  heat  flow  may  be 
applied.   The  following  development  has  been  principally 
adapted  from  Harrach  [30,32],  but  other  authors  have  similar 
developments  [7,49]. 

Suppose  that  at  time  t  =  0  laser  radiation  is  directed 
upon  an  opaque,  solid  slab  of  metal.   A  portion  of  the  radia- 
tion, is  absorbed;  the  rest  being  reflected  as  described  in 
the  previous  section.   A  temperature  distribution  T(r,t) 
develops  throughout  the  material.   When  a  temperature  gradient 
exists  in  a  body,  there  is  an  energy  transfer  from  the  high- 
temperature  region  to  the  low- temperature  region.   The  heat 
conduction  equation  (without  phase  changes)  is  [15] 


V-J(r,t)  +  PC  9T^,t)   =   A(r,t)  (1) 
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where  T(r,t)  is  the  temperature,  J(r,t)  is  the  thermal  energy 
crossing  unit  area  per  unit  time,  pc  is  the  heat  capacity  per 
unit  volume,  and  A(r,t)  is  the  net  heat  energy  per  unit 
volume  time  generated.   Fourier's  law  relates  the  temperature 
to  the  heat  flux 

J(r,t)   =   -KVT(r,t)  (2) 

To  simplify  the  calculations,  it  is  assumed,  as  in   Refs.  4  and  59, 
that  the  problem  is  planar,  heat  is  introduced  only  from  the 
surface  and  flows  into  the  material  in  the  x-direction.   The 
thermophysical  constants  are  independent  of  temperature, 
radiation  and  convection  from  the  surface  is  negligible,  and 
the  liquid  phase  of  the  material  can  be  ignored.   The  heat 
flow  can  now  be  modeled  using  the  Fourier  heat-conduction 
equation  [54] 

32T(x,t)    1   ,       =   1  JT(x,t)  , 

2    +   Aix,t;      k    5  u; 

3x 

where  k  is  the  thermal  diffusion  constant  equal  to  pKC  and 
K  is  the  thermal  conductivity  constant.   The  rate  of  heat 
production  can  be  written  as 


A  =   —  (4) 

A      pC6  [q) 


where  a  is  the  surface  absorptivity,  I  is  the  flux  density 
of  the  incident  radiation,  5  is  the  penetration  skin  depth  of 
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radiation  into  the  target  material,  p  is  the  density  of  the 
target  material,  and  C  is  the  specific  heat.   Now  by  substi- 
tution of  (4)  into  (3),  the  model  equation  becomes 


.  92T(x,t)    3T(x,t)      ctl(x,t) 

K   I  2   +   at        pcs 

3x 


In  order  to  solve  this  differential  equation,  the  various 
parameters  must  be  known  or  at  least  estimated.   With  the 
exception  of  6,    each  parameter  has  an  accepted  value  which 
can  be  used.   The  value  for  5  is  not  so  easy  to  determine. 
It  is  expected  to  be  small,  10    cm  [74].   It  can  be  assumed 
that  the  laser  radiation  is  absorbed  in  an  inf initesimally 
thin  surface  layer  (5  ■*  0)  .   The  source  term  may  now  be 
deleted  from  (5) .   The  term  is  not  lost  however,  since  it 
will  be  incorporated  into  the  surface  boundary  conditions. 
Equation  (5)  now  becomes 


_  92T(x,t)    3T(x,t)   _   -  ,,. 

-  K  ——2 +  — ^ -   0  (6) 

3x 


The  specific  boundary  conditions  depend  on  whether  the  sur- 
face of  the  target  is  being  vaporized  or  not.  For  the  case 
at  hand,  no  vaporization,  the  surface  boundary  condition  is 


„  JT  x,t)      al  t    -        _ 
-  K  — tt-5 — - — -  =   — —-^   for  x  =  0 

9x  pC 
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or 


K  dl^X,t)       =      al(t)  (7) 


A  second  boundary  condition  arises  when  the  target  is  con- 
sidered to  be  a  semi-infinite  solid. 


lim  3T(x,t)   _   n  ,Q, 


And  finally,  the  third  condition  is 

T(x,0)   =   0  (9) 

In  general,  equation  (6)  does  not  yield  an  exact  solution. 
One  method  of  solution,  utilized  by  Harrach,  is  the  "heat 
balance  integral  method"  [30].   This  method,  which  is  based 
on  the  assumption  that  "the  partial  differential  equation  is 
required  to  be  satisfied  only  in  an  average  sense  throughout 
the  solid  rather  than  at  each  point  x,"  is  outlined  below. 
The  integral  method  reduces  the  nonlinear  boundary  value 
problem  to  an  ordinary  initial  value  problem  where  solutions 
can  be  expressed  in  closed  analytical  form. 

Integrating  eqution  (6)  over  the  spatial  interval  of  the 
target,  x  _<  x  <_  I   gives 
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2  1 

r   k  — 2  dx   =     /  at  ^ 

x  '     9x         x 
s  s 


In  general  taking  the  derivative  of  an  integral  yields 


ft    /   TdX   "     /f  dx  +  T(x-l.t)  || 


s  s 


di: 

dx 


-  T(x=  x  ,t)  -r^- 
s     at 


Using  equation  (10)  and  the  results  of  equation  (11) ,  noting 
d£/dt  =  0  and  assuming  k  is  constant 


k(^) 


X=£ 


_T 
x 


,      £  dx 

)   -|-    /  T  dx  +  T(xs,t)  ^     (12: 

X=X  X 

s  s 


If  the  target  surface  is  stationary  then  dx  /dt  =  0.   Goodman 
[26]  points  out  that  the  key  step  in  the  careful  choice  of 
the  solution  form,  T(x,t).   The  idea  is  to  choose  T(x,t)  such 
that  it  can  be  integrated  explicitly.   This  gives  an  ordinary 
first  order  differential  equation  when  substituted  into 
equation  (12) .   Now  to  apply  this  technique  to  the  problem 
defined  by  equation  (6)  and  the  boundary  conditions  (7) , 
(8) ,  and  (9) 


Sr  /  T(x,t)dx  =  JSlitl.  (13) 

dt  0'  K 
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Consider  the  solution  form 


T(x,t)   =   T(0,t)  [1  -x/f(t)]2  exp[-x/c(t)]  ,  (14) 


for  0  <  x  <  E,  (t) 


and 


T(x,t)   =   0   for  x  >  £  (t) 

where  T(0,t)  is  the  front  surface  temperature  and  £(t)  is  a 
time  dependent  thermal  penetration  depth.   Integrating 
equation  (14)  over  x  from  0  to  £  and  incorporating  equation 
(13)  yields  the  ordinary  differential  equation 


g£[T(0,t)5(t)]   =   [kaI(t)]/[l-2e_1]K  (15! 


Substituting  equation  (14)  into  equation  (7)  gives 


5(t)   =   3KT(0,t)/aI(t)  (16) 


Solving  equations  (15)  and  (16)  simultaneously  yields 


T(0,t)   =   ( ^ ~)1/2  (al(t)   /  I(t')dt')  1/2 

3(l-e2  X)K  0 

(17) 
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and 


5(t)   =   (_3k   ,1/2  ,  ;*  cl(t-)df  ,1/2         (18) 
l-2e        0      aj-^ 


Several  special  cases  can  now  be  considered  by  stipu- 
lating the  laser  profile  I(t)  and  the  dependence  of  the 
absorptivity  on  temperature,  a  =  a(t). 

CASE  1: 

Let  the  absorptivity  coefficient  be  constant,  a  ,  and 
the  laser  pulse  be  a  step  function.   The  solution  to  equation 
(17)  is  now  found  to  be 


T(0,t)   =   ^( ^-y-)1/2  (19) 

3(l-2e  x) 


and  from  equation  (18) 


C(t)   =   (3KT/(l-2e_1)1/2  (20] 


For  comparison,  the  exact  solution  as  outlined  by  Carslaw 
and  Jaeger  [61]  is 


T(x,t)      =     — ^(— )1/2{ ^t-tk  exp(-x2/4kt) 

K        *  4Ukt)1/2 


(21) 


-  |erfc[x/(4^t)1/2] } 


In  the   above   equation   erfc  denots   the   complementary   function 
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erfc(X)   =   1  -  erc(X) 


v 


/       2 
=   2/tt   /  exp(-y  )dy 


The  agreement  between  the  approximation,  equation  (19)  and 

the  exact  solution,  equation  (21)  seems  to  be  quite  good. 

For  T(0,t),  equation  (21)  is  less  than  1%  higher. 

Now  the  solution  just  described  applies  up  to  a  time  t 

At  T(o,t  )  =  T   the  vaporization  time  is 
v     v       r 


t    =   {3(1  -  2e-1)/kHKT  /a    I  }2  (22) 

V  V   o  o 


for  I  /3KT   >  1. 
o    v  — 

CASE  2: 

The  laser  pulse  is  a  step  function  and  the  absorption 
coefficient  increases  linearly  with  the  surface  temperature 
up  to  the  vaporization  temperature  T   according  to  the 
relation 


a   =   [a+  (a  -  aQ)  }   [T(0,t)/Tv] 


The  solution  of  the  heat  flow  equation  (14)  at  x  =  o   is 
then  given  by 


T(0,t)   =   (a  -  ao)Tv/av-dQ  (23) 
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with 


?(t)   =   3KTQ(a  -  dQ)/[a(ctv-  ao)IQ]  (24) 


and  t   is  as  before, 
v 

C.   MATERIAL  REMOVAL 

Material  removal  begins  with  vaporization.   As  vaporiza- 
tion proceeds,  the  absorption  of  the  laser  beam  increases 
rapidly  as  the  depth  of  the  crater  increases.   At  first  the 
evaporation  rate  and  boundary  temperature  are  equal  to  zero 
and  all  of  the  laser  radiation  is  used  in  heating  the  opaque 
solid  as  pointed  out  in  the  previous  section.   If  the  external 
flux  is  intense  enough,  the  surface  will  be  brought  to  the 
vaporization  temperature  T   at  time  t  .   For  later  times, 
t  >  t  ,  the  problem  is  still  defined  by  equation  (6) ;  but 
one  of  heat  transfer  in  the  presence  of  a  moving  phase  boun- 
dary requiring  new  boundary  conditions. 

For  the  first  boundary  condition,  the  surface  boundary 
energy  balance  relation,  equation  (7)  must  be  changed  to 
include  a  term  to  account  for  the  energy  expended  on  melting 
and  vaporization.   It  must  also  account  for  the  fact  that 
the  surface  is  no  longer  stationary. 


k  3T(x,t) 


3x 


dx 

=   txl(t)  -  pl  ^-  (25) 

/  j.\  Qt 

x=x  (t) 
s 
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where  x  (t)  is  the  instantaneous  position  of  the  surface. 
s 

For  example,  x   =  0  when  t  <  t  .   L  is  the  effective  latent 
c  s  —  v 

heat  of  sublimation  described  by  L  =  L   , .  +  L      .   The 

2  melt     vapor 

second  boundary  condition  is  as  before,  that  is  the  target 
is  assumed  to  be  a  semi-infinite  solid 


lim  3T(x,t)   =   o  (26) 

a  X 

x->°° 


The  final  condition  is  based  on  the  temperature  distribution 

in  the  material  at  t  . 

v 

T(x,t)   =   f(x)  (27) 

where  the  function  f(x)  is  determined  by  the  prevaporization 
solution. 

With  the  boundary  conditions  just  described,  equation 
(6)  becomes 


j      *  i  j  /±.\  dx  (t) 

k<  /  TdX,  .  |ui(t,  -  PL  Sfja.]   -  Tv  -fc- 

x 

s 


If  dx  /dt  =  0  then  equation  (28)  reduces  to  the  prevaporiza- 
tion  form,  equation  (13) ,  as  expected.   The  trial  solution 
for  the  vaporization  case  may  be  chosen  to  be 


x-x  (t)     2       x-x  (t) 
T(x,t)   =   Tv[l  -  s(t^.x  (t)]   exp[-g(t)fx  (t))] 


(29) 
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for  0  <  x-x  (t)  <_  5(t)  -x  (t)  and  T(x,t)  =  0  for 

x-x  (t)  >  ?  (t)  -x  (t)  .   The  surface  temperature  is  fixed 

at  T   as  the  position  of  the  surface  x  (t)  moves  into  the 
v        r  S 

solid.   The  surface  recession  velocity  can  be  represented 
by 


ct(T  )I(T  )    .  ,     CT 

x  (t)   =   ^ —t(ir-)  ~  1]{  [1-4(3  +  -T1)]   [(—)  -  1]  } 

S  pL         t  z       L       t 

K  V  V 

(30) 
When  the  material  is  exposed  to  large  constant  flux  and 
begins  vaporizing  after  time  t  ,  the  rate  of  material  removed 
will  approach  a  steady  state  as  the  corresponding  surface 
recession  velocity  reaches  a  steady  state. 


X     =   a(T  )I(T  )/p  (L  +  CT  )  (31) 

SS  V '   v  v  V 


The  steady  state  penetration  depth,  D(t)  can  now  be  written 
as 


D     =   kpL/KT    =   L/CT 
ss  v         v 


Figure  7  from  Ref.  30  depicts  the  normalized  thermal  penetra- 
tion depth  D(t) ,  surface  position  x  (t) ,  and  surface  recession 
velocity  x  (t)  for  aluminum.   Next  the  crater  depth  l      at 
a  specific  time  t   can  be  determined  from  the  integral 
equation 
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Figure  7.   Time  dependence  of  normalized  values  of  surface 
recession  velocity  x  (t) ,  thermal  penetration 
depth  D(t),  and  surface  position  x  (t) ,  for 
aluminum  from  Ref.  30. 
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£         /    x  (t)  dt  (32) 

c      .  '  s 
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assuming  that  the  target  is  effectively  infinitely  thick. 

The  previous  discussion  indicates  how  the  temperature  of 
a  local  spot  on  a  target  can  be  raised  to  the  point  where 
vaporization  begins.   At  the  same  time,  thermal  diffusion 
begins  as  the  local  spot  temperature  rises.   As  a  result, 
for  most  metals,  the  material  is  removed  as  a  consequence 
of  three  processes:   vaporization,  the  development  of  pressure 
in  the  cavity  by  vapor  expansion  and  the  creation  of  molten 
metal.   Generally,  vaporization  begins  at  the  surface.   The 
initial  velocity  of  a  vapor  jet  flowing  into  a  vacuum  is 
equal  to  the  local  sound  velocity.   As  the  crater  depth 
increases,  the  vapor  velocity  increases  to  a  supersonic  flow 
[4],  which  provides  a  mechanism  for  the  washing  of  liquid 
metal  from  the  walls  of  the  crater.   The  fraction  of  material 
removed  in  the  liquid  state  increases  with  pulse  duration  [21] 
because  of  the  increase  in  temperature  of  the  crater's  inner 
walls  with  time. 

Theoretical  results  are  generally  difficult  to  compare 
to  experimental  ones  because  of  several  factors.   First,  the 
amount  of  molten  metal  washed  out  of  the  crater  and  the 
initial  fraction  of  absorbed  energy  is  not  well  known.   Next, 
a  change  in  absorbed  energy  with  temperature  is  generally 
unknown.   And  finally,  estimates  of  the  area  on  which  the 
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beam  was  focused  can  involve  considerable  error.   The  dis- 
tribution of  illumination  over  a  spot  can  be  very  nonuniform. 

In  addition,  experimental  results  on  the  amount  of  material 
vaporized  by  a  normal  pulse  laser  exhibit  considerable  varia- 
tion even  for  reported  "similar  conditions".   The  literature 
provides  a  number  of  sources  where  values  of  the  amount  of 
metal  ejected  by  a  normal  pulse  laser  are  given  [59,58,11, 
28,4,21,2] . 

Figure  8  depicts  the  ejected  mass  as  a  function  of  beam 
energy  for  a  700  us  duration  pulse  from  a  normal  pulse  neo- 

dymium  glass  laser  [11].   The  beam's  total  energy  output  was 

-3    2 
1-15  joules  and  its  spot  size  was  10    cm  .   In  each  case, 

material  removal  was  accompanied  by  breakdown.   The  group 

of  low  melting  metals  experience  larger  amounts  of  material 

removed.   Figure  8  shows  that  the  amount  of  material  ejected 

increases  rapidly  (exponentially)  with  increased  energy. 

Chun  and  Rose  conducted  similar  experiments  a  few  years 

later  [21] .   Figure  9  shows  the  dependence  of  material  removal 

on  pulse  duration  at  an  average  laser  power  of  30  KW.   To 

obtain  the  data  in  Figure  9,  they  used  different  laser  pulse 

durations  but  constant  average  laser  power.   The  mass  of 

the  ejected  material  was  determined  by  measuring  the  amount 

of  material  collected  on  the  walls  of  a  glass  box  surrounding 

the  specimen. 

In  conclusion,  one  can  say  that  the  effect  of  normal 

pulse  laser  radiation  on  the  surface  of  an  opaque  solid  is  a 
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Figure  8.   Experimental  mass  removal  for  various  metals 
using  a  700  ms  Nd-glass  laser  focused  to 
10"3  cm2  (Ref .  11) . 
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Figure  9.   Time  dependent  material  removal  measurements. 
These  measurements  are  the  average  weight  of 
material  removed  from  three  specimens  for  _     _ 
laser  pulses  having  an  average  power  of  10  w/cm 
but  different  time  durations  (Ref.  21). 
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pronounced  one.   I  have  considered  an  ordinary  thermodynamic 
approach  for  the  calculation  of  the  target  temperature  dis- 
tribution, time  to  vaporization,  and  crater  depth.   Experi- 
mental results  on  material  removal  for  lasers  providing 
fluxes  up  to  10   -  10    watts  have  been  outlined.   Repre- 
sentative values  for  mass  removed  and  crater  depth  are  on 
the  order  of  a  few  milligrams  and  a  few  millimeters, 
respectively . 
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III.   EXPERIMENTAL  DESIGN 

A.   EQUIPMENT 

The  arrangement  of  the  equipment  utilized  during  this 
experiment  is  depicted  in  Figure  10. 

1.  Laser  System 

A  neodymium  glass  laser  emitting  a  wavelength  of  1.06 
micrometers  was  used  for  this  experiment.   The  pockel ' s  cell 
of  the  Korad  K-1500  laser  was  replaced  by  a  simple  mirror  to 
allow  for  normal  pulse  operation.   Figure  11a  shows  a  typical 
oscilloscope  display  of  the  laser  pulse.   Figure  lib  is  an 
expanded  view  of  the  same  pulse  to  show  the  individual 
spikes.   A  detailed  description  of  the  laser  system  can  be 
found  in  Appendix  A  of  Ref.  55. 

2.  Target  Chamber 

The  target  chamber  was  made  from  a  6  inch  cube  of 
unbaked  aluminum.   Its  internal  volume  is  12.9  +  .03  liters. 
The  vacuum  system  using  a  mechanical  forepump  and  an  oil 
diffusion  pump,  was  capable  of  providing  chamber  pressures 
on  the  order  of  10  "  Torr.   The  target  and  ionization  gauge 
were  mounted  as  shown  in  Figure  12.   The  targets  were  made 
of  6062  aluminum,  50  mm  in  diameter  and  7  mm  thick.   The 
surfaces  were  machined  "smooth"  on  a  metal  lathe.   A  smaller 
target  made  from  sheet  aluminum  (12mm  x  20mm  x  .5mm)  was  also 
used.   They  were  mounted  on  a  holder  described  by  Hwang 
[35]  . 
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Figure  10.   Block  diagram  of  the  general  experimental 
equipment  layout. 
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Figure  11a. 


Oscilloscope  display  of  the  laser  pulse  at 
100  us/div. 


Figure  lib. 


Oscilloscope  display  of  the  laser  pulse  showing 
the  individual  spikes.   Time  is  2  ys/div. 
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Figure  12. 


Top  view  of  the  target  chamber.   There  is  a 
port  in  the  top  of  the  chamber  where  the 
detector  for  the  plasma  radiation  was 
located.   C  denotes  the  positive  ion 
collector,  G  denotes  the  accelerating  grid 
and  F  the  Filament. 
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3 .  Ionization  Gauge 

A  Hughes  ionization  gauge  tube  type  6578  was  used  as 
a  device  to  measure  the  material  ejected  from  the  surface  of 
the  target.   Gauge  #1  was  placed  inside  the  target  chamber 
in  the  path  of  the  excited  material.   Gauge  #2  was  located 
below  the  target  in  the  manifold  to  the  chamber.   Figure  13 
pictures  the  gauge  #1  plus  mounting  hardware.   The  gauge  was 
designed  with  a  very  fine  collector  wire  suspended  in  a 
cylindrical  grid.   A  filament  placed  outside  the  grid  is 
heated  by  passing  a  current  (6  amps)  through  it.   The  grid 
cylinder,  which  has  a  positive  potential  applied  to  it, 
accelerates  the  electrons  from  the  filament.   A  simple  elec- 
tric circuit  for  an  ionization  gauge  is  diagramed  in  Figure  14. 
There  is  a  high  probability  that  collisions  between  electrons 
and  molecules  will  result  in  the  formation  of  positive  ions. 
Since  the  collector  has  a  negative  applied  voltage,  it  collects 
ions  within  the  grid  area.   The  variations  of  collector  cur- 
rent can  then  be  measured  and  equated  to  the  number  of  par- 
ticles in  the  chamber  [23]  .   If  a  plasma  passes  through  the 
ion  gauge,  the  collector  will  sense  the  positive  ions  resulting 
in  an  appropriate  current  increase.   This  is  similar  in  opera- 
tion to  an  electric  probe  used  for  plasma  diagnostics  [33]. 

4 .  Instrumentation 

The  laser  output  signal  was  sensed  using  a  Lasermetric 
Series  3117,  High  Speed  Detector.   To  protect  the  detector, 
50,  5,  and  1  percent  transmission  filters  (at  1.06  ym)  were 
placed  in  line  with  the  radiation.   The  signal  from  the 
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Figure  13. 


Hughes  Ionization  Gauge,  gauge  1,  with 
mounting  hardware. 
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Figure  14.   An  example  of  a  sample  electrical  circuit 
for  an  ionization  gauge.   C  denotes  the 
positive  ion  collector,  G  denotes  the 
accelerating  grid,  and  F  the  filament. 
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detector's  high  speed  (<  150  picosecond)  PIN  diode  was  dis- 
played on  a  Tektronix  7904  Oscilloscope.   Total  energy  in 
the  laser  pulse  was  determined  using  a  Laser  Precision  RK- 
3200  Series  Pyroelectric  Energy  Meter.   A  10  percent  trans- 
mission filter  (at  1.06  pm)  was  placed  in  front  of  the  detec- 
tor face,  to  insure  that  the  level  of  intensity  reaching  the 
detector  did  not  exceed  the  damage  threshold  of  the  detector. 
This  instrument,  which  gave  direct  digital  display  of  the 
total  energy  in  a  pulse,  was  calibrated  using  a  Hadron  Model 
118  Thermopile.   The  calibration  accuracy  was  no  better  than 
20  percent  which  is  to  be  expected  [59]. 

The  plasma  radiation  was  detected  using  a  Hewlitt 
Packard  pin  diode  detector  with  response  time  <  1  nanosecond. 
To  eliminate  any  signal  from  1.06  ym  radiation,  a  filter  with 
optical  density  of  4.5  at  1.06  um  was  placed  over  the  face 
of  the  detector.   The  detector's  signal  was  displayed  on  a 
Tektronix  7704  Oscilloscope. 

The  two  ion  gauges  in  the  system,  one  in  the  chamber 
and  one  in  the  elbow  beneath  the  chamber,  were  each  connected 
to  a  Veeco  Vacuum  Gauge  Control  Panel.   An  oscilloscope  was, 
in  turn,  connected  to  the  recorder  connection  on  each  panel. 
The  amplifier  output  was  -2  volts  with  respect  to  the  ground 
for  a  full  scale  pressure  reading  on  all  scales.   A  set  of 
Granville-Phillips  control  panels  was  also  used  in  lieu  of 
the  Veeco  models.   These  control  panels  provided  a  -100 
millivolt  output  for  a  full  scale  pressure  reading  on  all 
scales. 
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Due  to  the  singularity  and  duration  of  the  pulses, 
all  oscilloscope  displays  had  to  be  recorded  by  photographing 
with  the  appropriate  Tektronix  camera.   Polaroid  3000  speed 
black  and  white  film  was  used  exclusively. 

Target  weight  measurements  were  made  with  a  Cahn 
Electrobalance,  Model  DTL.   This  scale  was  capable  of 
measurements  as  small  as  .01  milligrams  with  a  reported 
accuracy  to  the  nearest  .01  milligram. 

B.   PROCEDURE 

1.  Ion  Gauge  Measurement  Study 

A  target  was  affixed  to  the  target  holder  and  the 
chamber  evacuated  to  about  10  "  Torr.   Prior  to  taking  data, 
all  oscilloscopes  and  ionization  gauges  were  adjusted  and 
calibrated. 

For  all  experiments  the  laser  was  set  to  provide  an 
output  of  15-18  Joules.   The  beam  intensity  (focus)  was  ad- 
justed with  a  simple  lens  focusing  system.   A  complete  sequence 
of  shots  was  taken  varying  the  influence  on  the  target,  while 
measurements  of  the  plasma  and  laser  outputs  and  ionization 
gauge  outputs  were  made.   Some  shots  were  repeated  with 
gauge  filament  both  off  and  on  to  distinguish  the  effect  of 
the  plasma  on  the  collector  current  and  to  show  the  effect  of 
the  ionizing  filament  electrons. 

2.  Material  Mass  Study 

The  system  was  prepared  and  operated  as  before  except 
that  special  small  targets  were  used.   The  procedure  began 
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by  weighing  a  "standard"  on  the  Electrobalance  and  recording 

the  weight.   The  target  sample  was  then  weighed.   After  5 

or  6  laser  shots  on  the  target  with  the  beam  focused,  it 

was  reweighed.   The  "standard"  was  also  reweighed  to  determine 

if  the  balance  had  "drifted".   The  difference  in  weight 

readings  for  the  sample  as  well  as  the  "standard"  were 

recorded. 
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IV.   RESULTS  AND  DISCUSSION 

One  objective  of  this  research  was  to  find  a  method  by 
which  one  could  measure  the  amount  of  material  ejected  from 
a  target  irradiated  by  a  laser  beam.   Initially,  a  Granville- 
Phillips  ionization  gauge  control  panel  was  used  to  monitor 
the  ionization  gauge  signals,  but  was  found  to  be  unacceptable 
The  control  panel's  electronic  circuits  saturated  when  the 
laser  was  fired,  yielding  outputs  that  were  not  characteristic 
of  the  chamber  pressure  fluctuations.   The  data  described 
was  collected  using  Veeco  control  panels.   Through  the  use 
of  an  ionization  gauge  installed  in  the  target  chamber, 
gauge  1,  I  was  able  to  sense  the  plasma,  neutral  particles, 
and  desorbed  gases  evaporated  from  the  target  surface.   The 
reaction  of  the  ionization  gauge  to  each  will  be  discussed 
below. 

A.   PLASMA 

When  the  600  ps  pulse  was  focused  on  an  aluminum  target, 
material  was  evaporated  and  a  crater  formed  on  the  surface. 
The  crater  was  oblong  having  dimensions  of  2mm  by  1mm  by 
0.5mm  deep  for  pulses  with  energy  15-18  J.   Assuming  a  coni- 
cal shape  for  the  crater,  the  calculated  mass  removed  is 
0.71  mg;  however,  not  all  of  this  material  is  vaporized. 
The  crater  was  ringed  with  metal  ejected  from  the  crater. 
This  metal  was  deposited  on  outer  edges  as  a  result  of  the 
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vapor  flow  scraping  molten  metal  off  of  the  inner  cavity 
walls.   A  detailed  display  of  the  laser  output  and  target 
plasma  radiation  from  the  same  pulse  is  depicted  in  Figure  15. 

The  beam  was  focused  resulting  in  an  energy  density  approxi- 

3      2 

mately  1.9  xlO   J/ cm  .   Light  emission  from  the  plasma  was 

measured  with  a  photo  detector  located  on  the  top  of  the 
target  chamber.   A  neodymium  absorption  filter  was  used  to 
eliminate  any  direct  radiation  from  the  laser.   Figure  15 
indicates  that,  for  this  case,  all  of  the  laser  spikes  result 
in  a  plasma  being  generated.   The  laser  spike  is  relatively 
symmetrical,  but  the  plasma  emission  lineshape  rises  sharply 
as  breakdown  threshold  is  reached  and  then  trails  off  exponen- 
tially as  the  plasma  expands,  cools,  and  recombines.   There 
is  also  a  good  correlation  between  plasma  radiation  and  the 
ionization  gauge  output.   Figure  16  depicts  these  two  traces 
for  a  single  laser  pulse  which  was  defocused  more  than  the 
pulse  in  Figure  15.   Figure  16a  shows  breakdown  spikes  at 
t,  =  0.38  ms  and  another  at  t_  =  0.66  ms  with  approximately 
the  same  amplitude.   The  corresponding  ionization  gauge  output 
can  be  seen  in  Figure  16b.   The  two  plasma  spikes  at  t,  are 
manifested  as  one  ionization  gauge  deflection.   Both  curves 
show  an  exponential  decay.   Of  particular  interest  is  the 
0.8  us  amd  24  ys  time  for  the  1/e  decay  of  the  plasma  radia- 
tion and  ionization  gauge,  respectively.   When  the  gauge's 
filament  was  turned  off,  the  signal  displayed  on  the  oscillo- 
scope was  the  same.   This  is  to  be  expected  because  the 
negatively  biased  collector  acts  like  an  ion  collecting 
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a.   Oscilloscope  trace  of  laser  output  at  1.06  ym, 


Oscilloscope  trace  of  plasma  radiation  with  1.06  ym 
narrowband  absorption  filter  in  place. 


Figure  15.   Oscilloscope  traces  of  laser  output  and  plasma 
radiation  taken  from  the  same  pulse.   Trigger 
was  delayed  0.22  ms .   The  pulse  energy  was  15 
Joules.   The  beam  was  focused  to  approximately 
8  x  10~3  cm^ .   Time  reference  for  both  traces 
was  5  ys/div. 
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Oscilloscope  trace  of  plasma  light.   The  wide  base 
signal  is  due  to  laser  flash  lamp  radiation.   The 
spikes  above  the  background  noise  are  light  emitted 
from  a  breakdown  plasma. 
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b.   Oscilloscope  trace  of  chamber  ionization  gauge  output  with 
filament  on  and  control  panel  set  on  the  10  '  Torr  range. 

Figure  16.  Oscilloscope  traces  of  plasma  radiation  and 
ionization  gauge  output.  The  time  for  both 
traces  was  0.1  ms/div. 
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electric  probe.   The  plasma  contains  many  ions  which,  in 
turn,  provide  a  signal  for  the  ionization  gauge.   Figure  16 
also  shows  that  there  was  no  detectable  time  lag  between 
plasma  radiation  and  the  corresponding  ionization  gauge 
signal.   Actually,  the  plasma,  traveling  at  about  2  x  10   cm/s 
[14],  took  approximately  5  us  to  cross  the  10.4  cm  distance 
from  the  target  to  the  collector,  but  this  time  was  not  per- 
ceivable on  the  0.1  ms  scale  in  Figure  16. 

B.   DESORBED  GASES 

To  demonstrate  that  the  ionization  gauge  deflections 

measured  "fly  by"  material,  the  laser  beam  was  defocused  to 

-2    2 

approximately  1.8  x  10    cm  .   At  this  point,  there  was 

2 
insufficient  energy  density,  J/cm  ,  to  ionize,  melt,  or 

vaporize  the  target  surface.   Only  gases  were  desorbed  off 
of  the  surface.   A  visual  inspection  of  the  target  showed 
a  0.4  cm  diameter  circle  where  the  target  had  been  heated. 
Figure  17  shows  traces  from  the  two  different  ionization 
gauges  mounted  in  the  system.   The  upper  trace  is  from  the 
chamber  gauge,  gauge  1,  while  the  lower  one  corresponds  to 
the  gauge  located  below  the  chamber,  gauge  2.   The  top  trace 
shows  a  definite  fluctuation,  while  the  lower  trace  is 
relatively  flat  though  its  gain  was  5  times  greater.   When 
the  laser  is  fired  several  times  in  the  same  spot,  the  first 
three  or  four  pulses  caused  a  much  higher  positive  pressure 
indication  than  subsequent  ones.   This  effect  has  been  pre- 
viously reported  [55,62].   It  is  due  mainly  to  the  desorption 
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Figure  17.   Oscilloscope  display  of  the  ionization  gauge 
output  with  the  control  panel  range  set  at 
10~"  Torr  range.   The  upper  and  lower  displays 
are  from  gauge  1  and  2  respectively.   The  beam 
was  defocued  to  1.8  x  10~2  cm2  (no  plasma 
generation).   Time  was  0.2  ms/div.   Gain  was 
1  V/div  for  the  upper  trace  and  0.2  V/div. 
for  the  lower. 


Figure  18. 


Oscilloscope  trace  of  #1  ionization  gauge.   A 
700  G  horseshoe  magnet  was  placed  parallel  to 
the  target  surface  to  deflect  charged  particles 
form  the  target  surface.   The  beam  with  energy 
of  15  J  was  focused  to  approximately  1.8  x  10~2 
The  reference  for  time  was  0.2  ms/div.   Tn§6 
control  panel  was  set  to  operate  on  the  10 
Torr  range. 
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of  impurities  such  as  nitrogen,  oxygen,  hydrogen,  water,  and 
carbon  dioxide  found  on  the  target  surface. 

A  baffling  result  was  the  initial  signal  in  the  negative 
direction  with  the  onset  of  lasing.   If  a  magnet  is  placed 
adjacent  to  the  target  to  deflect  charges  particles,  the  nega- 
tive signal  disappears  and  only  a  positive  signal  is  seen, 
Figure  18.   If  the  gauge's  filament  is  turned  off,  the  trace 
is  flat  for  the  same  operaing  conditions.   No  voltage  drop 
on  the  +150  volt  accelerating  grid  was  observed.   One  explana- 
tion for  the  negative  signal  may  be  thermal  or  photo  electron 
emission  from  the  target  surface.   The  gauge's  filament  has  a 
negative  space  charge  around  it.   The  additional,  probably 
more  energetic,  electrons  from  the  target  surface  would 
increase  the  space  charge  and  act  as  a  shield  returning  (repel) 
thermal  electrons  back  to  the  filament  surface.   The  net 
effect  is  a  reduction  of  the  number  of  electrons  accelerated 
by  the  grid.   This,  in  turn,  means  less  ionization  takes 
place  and  therefore  a  drop  in  collector  current.   The  positive 
portion  of  the  trace,  starting  at  the  dips  minimum,  Figure  17, 
was  due  to  the  pressure  increase  by  gases  desorbed  from  the 
target  surface.   The  conclusion  that  can  be  drawn  is  that  the 
chamber  ionization  gauge  was  actually  sensing  what  was  flying 
past  it. 

C.   NEUTRAL  PARTICLES 

Neutral  particles  present  a  special  problem  of  measure- 
ment.  It  was  expected  that  electrons  from  the  ionization 
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gauge's  filament  would  ionize  a  fraction  of  the  neutral 
aluminum  particles  allowing  for  their  collection.   This,  in 
fact,  happened.   A  comparison  of  Figure  19  and  20  shows  that 
the  neutral  particles  (and  desorbed  gases)  were  not  collected 
if  the  filament  was  off,  resulting  in  a  steep  decline  and 
a  flat  trace  after  the  plasma  passed,  Figure  20. 

It  was  observed  that  when  breakdown  at  the  target  surface 

occurred,  plasma  and  neutral  particles  were  ejected.   The 

1/2 
plasma  velocity  is  on  the  order  of  (kT  /M . )    ,  approximately 

2  x  10   -  2  x  10   cms,  depending  on  the  laser  pulse  parameters 

[14] .   The  plasma  was  therefore  much  faster  than  the  neutral 

particles  which  leave  the  surface  at  the  thermal  velocity  of 

4 
2  x  10   cms  [59].   In  Figure  19,  this  is  demonstrated.   The 

plasma  passed  the  chamber  ionization  gauge  first,  immediately 
followed  by  the  neutrals.  The  velocity  of  the  neutral  parti- 
cles can  be  determined  from  the  times  depicted  in  Figure  21. 

2 
The  laser  beam  was  focused  to  approximately  0.018  cm  and  a 

700  G  horseshoe  magnet  placed  next  to  the  target  to  suppress 
thermal  and  photo  electrons  as  previously  discussed.   A  posi- 
tive rise  in  the  oscilloscope  trace  began  at  t~  =  0.42  ms , 
0.22  ms  after  lasing  began,  t-,.   The  first  arriving  neutral 

particles  traveled  the  10.4  cm  from  the  target  to  the  collector 

4 
at  5.2  x  10   cm/s.   The  most  probable  thermal  velocity  of 

an  aluminum  particle  leaving  the  target  surface  at  the  vapori- 

1/2 
zation  temperature  of  2330  °K  can  be  calculated  from  (2Kt/M)  '  . 

4 
It  is  3.8  x  10   cm/s.   The  velocity  of  particles  evaporated 
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^^^^ P^^^^'v^vS .'^1  "leutral  '  s 

Figure  19.   Oscilloscope  display  of  the  #1  ionization 
gauge  output  with  filament  on  and  with  the 
control  panel  set  on  the  10"^  Torr  range. 
Beam  was  focused  down  to  2  x  10" 2  cm2 . 
Time  is  0.5  ms/div. 


Figure  20.   Ionization  gauge  oscilloscope  trace  for  the 
same  conditions  as  those  in  Figure  17  except 
that  the  filament  current  was  turned  off. 
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Figure  21. 


Oscilloscope  trace  of  #1  ionization  gauge. 
A  700  G  horseshoe  magnet  was  placed  parallel 
to  the  target  surface  to  deflect  charged 
particles  from  the  target  surface.   The  beam 
with  energy  of  15  J  was  focused  to  approximately 
1.8  x  10    cm^ .   The  reference  for  time  was 
0.2  ms/div.   The  control  panel  was  set  to 
operate  on  the  10~5  Torr  range.       v1 
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from  the  surface  can  probably  be  characterized  by  a  Max- 
wellian  distribution.   The  faster  experimental  velocity 
value  can  be  attributed  to  particles  with  velocities  in  the 
tail  of  the  distribution. 

Figure  22  depicts  the  results  of  a  six  shot  sequence 
fired  at  one  spot.   Again,  a  magnet  was  used  to  deflect 

photon  and  thermal  electrons  away  from  the  filament.   The 

-2    2 

spot  size  of  1.8  x  10    cm   and  time  scale  of  0.22  ms/div 

were  kept  constant  as  the  laser  energy  was  increased  from 
18.8  to  36.6  Joules  (18.8,  19.6,  20.5,  22.9,  24.6,  and  36.6 
Joules) .   The  traces  have  increased  amplitude  with  increased 
energy.   The  relative  height  of  that  portion  of  the  curve 
attributed  to  neutral  particle  emission  from  the  target  was 
very  dependent  on  laser  flux.   For  all  cases,  positive  deflec- 
tions due  to  neutral  particles  occurred  at  the  same  time, 
t,  =  0.42  ms.   The  final  shot  shows  an  initial  spike,  caused 
by  a  small  plasma  cloud,  followed  by  a  sharp  rise  due  to 
neutral  particles.   One  can  conclude  that  as  energy  was 
increased,  the  initial  arrival  time  of  neutral  particles  at 
the  collector  remains  constant.   Further,  the  rise  time 
decreased  and  the  amplitude  of  the  pressure  signal  increased 
as  the  beam  energy  density  was  increased. 

As  an  example,  Figure  23a  depicts  the  ionization  gauge 
signals  as  material,  ejected  from  the  target  surface,  passed 
by  it.   In  this  case,  the  experiment  was  conducted  without 
the  target  magnet.   Figure  23b  is  the  corresponding  trace 
of  the  plasma  radiation.   Both  traces  were  from  the  same 
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Figure  22. 


Oscilloscope  trace  of  chamber  ionization 
gauge.   A  700  G  magnet  was  placed  parallel  to 
the  target.   The  laser  beam  was  increased  with 
each  shot  as  the  spot  size,  1.8  x  10~2  cm2, 
and  the  time,  0.2  ms/div.  were  kept  constant. 
The  gauge  control  panel  was  set  on  the  10~5 
Torr  range. 


63 


Oscilloscope  display  of  ionization  gauge  output  with 
control  panel  set  on  10~5  Torr  range.   The  lower  trace 
is  from  gauge  #2  with  the  filament  off.   The  small 
ripple  was  due  to  plasma  radiation.   The  upper  trace 
is  from  the  chamber  gauge. 


b.   Oscilloscope  trace  of  plasma  radiation. 

Figure  23.   Oscilloscope  displays  of  ionization  gauge 
outputs  and  plasma  radiation  for  the  same 
shot.   The  beam  which  was  defocused  to 
approximately  2  x  10~2  cm^  had  energy  of 
18  J.   The  time  reference  for  both  displays 
was  0.2  ms/div. 
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laser  pulse  and  had  the  same  time  scale  rf  0.2  r.s.  div.   rhe 

pulse,  with  energy  of  18  J,  was  recused  re  apprcxir.arely 

-2    2 
2  x  10    cm  .   With  the  lasing  starting  =:  0.22  r.s,  Figure 

23b  shows  a  piasr.a  spike  at  C.24  r.s.   The  upper  arc  lewer 

ionization  gauge  traces  aisc  shew  the  plasma  signal,  point 

a.   When  the  piasr.a  expands  and  cecls,  i.e.,  afrer  :essa:::: 

of  the  plasma  spike,  the  predominant  signal  in  rhe  negative 

direction  was  from  the  chore  :r  ch.err.al  elerrrrrs  erirrei 

from  the  laser  reared  surface,  segrer.r  b-c.   :'e;tral  parriri 

which  have  a  slower  velccity  char,  rhe  electrons,  arrived  ar 

the  gauge  last  causing  an  up-swir.g  in  rhe  rrace  ::  rhe  ngrr 

of  point  c.   Frr  a  case  rf  r.cre  incense  laser  radiation, 

consider  Figure  24.   The  rear,  had  ahrur  rhe  sar.e  energy, 

19.5  J,  but  ir  had  been  recused  r:  prrvide  breaker wn  as 

depicted  in  Figure  2  4b.   Bern  era  res  had  a  time  srale  rf 

0.2  ms/div.  to  facilirare  rrrparisrn.   The  gauge  crare  shrv: 

a  "wild"  deflection  as  rhe  piasr.a  passed  through,  segment 

continued,  still  adding  energy  cr  rhe  hot  surf are.   The  seg- 
ment b-c  shews  a  gradual  dselire  in  rhe  number  rf  neutral 
particles  ejected  frrr.  rhe  s urf are.   C reversely ,  when  las m: 

stops,  vaporization  steps  and  rhere  was  a  rapid  fall  in  rhe 
curve,  segr.enc  c-d.      "his  was  due  cr  rhe  rapid  surfare  rrrl: 

of  the  tarcet . 
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Oscilloscope  display  of  #1  ionization  gauge  output, 
Control  panel  range  was  set  on  10~5  Torr  range. 


b.   Oscilloscope  trace  of  plasma  radiation. 

Figure  24.   Oscilloscope  displays  of  ionization  gauge 
output  and  plasma  radiation  for  the  same 
shot.   The  beam  was  focused  from  that  in 
Figure  22.   Energy  was  approximately  18  J. 
Time  reference  was  0.2  ms/div.  for  both 
displays. 
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D.   MASS  OF  MATERIAL  REMOVED 

In  order  to  equate  the  chamber  ionization  gauge  signal 
to  a  specific  number  of  particles  ejected  or  ultimately  the 
mass  of  the  material  removed,  one  must  first  be  able  to 
measure  the  mass  of  the  ejected  material.   Small  aluminum 
targets  were  weighed  and  places  in  the  target  chamber.   The 

laser,  with  energy  about  18  J,  was  fired  at  five  locations 

2 
on  each  target.   Beam  spot  size  was  approximately  0.016  cm  . 

The  target  craters  were  approximately  0.5  mm  deep  and  had  the 
usual  metal  splashed  around  the  periphery  of  the  crater. 
From  the  volume  of  the  crater,  the  estimated  mass  removed 
was  0.7  mg.   However,  upon  trying  to  reweigh  the  target, 
it  was  found  that  only  a  few  hundredths  of  a  milligram  of 
material  seem  to  have  been  removed.   Compared  to  previous 
experimental  work,  this  amount  seems  to  be  extremely  low 
[11,21],   One  explanation  may  be  the  difficulty  encountered 
in  weighting  small  amounts.   The  Electrobalance  used  to  weigh 
the  targets  tended  to  drift.   Fluctuations  varied  from  hun- 
dredths to  tenths  of  milligrams  for  the  same  sample.   A 
standard  weight  was  used  to  eliminate  some  of  the  drift, 
but  only  the  long  term  effects  could  be  eliminated.   As  a 
consequence,  the  mass  of  material  removed  was  not  accurately 
determined. 
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V.   CONCLUSIONS 

Laser  induced  evaporation  of  material  from  a  target 
surface  has  been  investigated  using  a  normal  pulse  laser. 
The  technique  of  using  an  ionization  gauge  to  sense  the 
ejected  materials  does  work.   The  gauge  essentially  "sees" 
what  is  flying  past  it.   When  breakdown  on  the  target  surface 
occurs,  a  crater  is  formed.   The  breakdown  is  accompanied  by 
generation  of  a  plasma  and  neutral  particles.   Breakdown 
can  be  avoided  by  defocusing  the  laser  to  a  larger  spot  on 
the  target  surface.   In  this  case,  only  gas  desorption  from 
the  surface  occurs  and  much  smaller  signals  are  observed 
with  the  ionization  gauge. 

A  good  correlation  between  incident  radiation,  plasma 
radiation,  and  the  ionization  gauge  fluctuations  was  observed. 
For  laser  energy  densities  sufficient  to  cause  many  break- 
downs on  the  target  surface,  the  ionization  gauge  reacts  with 
large  positive  signals  to  the  plasma  which  arrives  first, 
and  then  to  the  slower  velocity  neutral  particles.   From 

time-of- flight  measurements  the  first  arriving  neutral  parti- 

4 
cle  velocity  was  determined  to  be  5.2  x  10   cm/s  which  com- 
pares favorably  with  published  data.   For  less  intense  laser 
radiation,  surface  breakdown  did  not  occur  or  occurred  only 
weakly  a  few  times  during  a  normal  laser  pulse.   For  this 
case,  the  ionization  gauge  initially  gave  a  negative  signal 
in  conjunction  with  the  onset  of  lasing  (excluding  positive 
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plasma  caused  fluctuations) .   This  can  be  attributed  to  ther- 
mal and  photo  electrons  ejected  from  the  target  surface. 
Once  the  slower  velocity  neutral  particles  (and  desorbed 
gases)  arrive  at  the  gauge,  the  signal  increases  and  becomes 
positive.   The  signal  amplitude  is  commensurate  with  the 
laser  energy  density. 

Mass  removal  measurements  were  inconclusive.   Experimental 
measurements  of  crater  volume  indicated  that  .017  mg  of 
aluminum  were  removed;  however,  some  of  this  material  was 
deposited  in  a  ring  around  the  crater  opening.   Further 
experimental  work  needs  to  be  done  on  weighing  the  mass 
ejected  from  the  target  surface.   Once  an  acceptable  technique 
is  developed,  then  correlation  experiments  between  the  weight 
of  the  mass  removed  and  ionization  gauge  fluctuations  should 
be  conducted  to  "calibrate"  the  ionization  gauge  output. 
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